Generalized covariant gyrokinetic dynamics of magnetoplasmas^ 

C. Cremaschini a , M. Tessarotto b ' c , P. Nicolini 6 and A. Beklemishev d 

"■Department of Astronomy, University of Trieste, Trieste, 
Italy, b Department of Mathematics and Informatics, 
University of Trieste, Trieste, Italy, 
c Consortium of Magneto-fluid-dynamics, 
University of Trieste, Trieste, Italy, 
d Budker Institute of Nuclear Physics, Novosibirsk, Russia 
(Dated: June 28, 2008) 

Abstract 

A basic prerequisite for the investigation of relativistic astrophysical magnetoplasmas, occurring 
typically in the vicinity of massive stellar objects (black holes, neutron stars, active galactic nu- 
clei, etc.), is the accurate description of single-particle covariant dynamics, based on gyrokinetic 
theory (Beklemishev et al., 1999-2005). Provided radiation-reaction effects are negligible, this is 
usually based on the assumption that both the space-time metric and the EM fields (in partic- 
ular the magnetic field) are suitably prescribed and are considered independent of single-particle 
dynamics, while allowing for the possible presence of gravitational/EM perturbations driven by 
plasma collective interactions which may naturally arise in such systems. The purpose of this work 
is the formulation of a generalized gyrokinetic theory based on the synchronous variational prin- 
ciple recently pointed out (Tessarotto et al., 2007) which permits to satisfy exactly the physical 
realizability condition for the four- velocity. The theory here developed includes the treatment of 
nonlinear perturbations (gravitational and/or EM) characterized locally, i.e., in the rest frame of a 
test particle, by short wavelength and high frequency. Basic feature of the approach is to ensure the 
validity of the theory both for large and vanishing parallel electric field. It is shown that the correct 
treatment of EM perturbations occurring in the presence of an intense background magnetic field 
generally implies the appearance of appropriate four-velocity corrections, which are essential for 
the description of single-particle gyrokinetic dynamics. 

PACS numbers: 52.30.Gz,45.20.Jj,52.27.Ny 
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I. INTRODUCTION 



From the mathematical viewpoint, it is well known that the Hamiltonian dynamics 
of charged single particles can be expressed, in principle, in arbitrary hybrid (i.e., non- 
canonical) state variables. The search of phase-space transformations yielding "simpler" 
equations of motion has motivated in the past theoretical research in physics and math- 
ematical physics. Among such transformations, a particularly significant case is provided 
by the gyrokinetic transformation, yielding generally hybrid variables defined in such a way 
that one of them, the so-called gyrophase angle, results ignorable [l, ^]. This transfor- 
mation, which is typically based on a perturbative expansion (the so-called Larmor radius 
expansion), can only be introduced when particles move in a sufficiently strong magnetic 
field. The gyrokinetic description of particle dynamics is useful because in these variables 
the equations of motion for single charged particles generally exhibit a significantly lower 
computational complexity (with respect to the case of non-gyrokinetic variables). This is 
due to the fact that the dependence on the fast gyro-angle formally disappears in this rep- 
resentation. The same approach, nevertheless, can be adopted in principle also for the 
formulation of the kinetic and fluid descriptions of relativistic plasmas, provided the treat- 
ment of these plasmas can be based on single-particle dynamics, i.e., radiation effects (such 
as, in particular, radiation-reaction effects) are negligible. For a certain class of problems in 
plasma-physics and astrophysics, involving the study of relativistic plasma flows in strong 
gravitational fields, the adoption of the gyrokinetic formulation results of basic importance. 
This involves the description of plasmas in which the strength of the magnetic field becomes 
very large, which are observed or assumed to exist in accretion disks and related plasma 
jets around neutron stars, black holes, and active galactic nuclei 3), Q]. The conventional 
non-relativistic gyrokinetic theory (GKT) has been previously generalized to include only 
some relativistic effects, which are essential for description of confinement of fusion products 
in thermonuclear reactors, namely, the particle itself is considered relativistic, while its drift 
velocity is not jf|. Other approaches have instead investigated the formulation of covariant 
relativistic theories applicable in astrophysical problems. A first attempt at a formulation 
of this type was the Ph.D. thesis of Boghosian [5J. In line with previous non-covariant treat- 
ments, this approach assumed however a flat space-time, a vanishing parallel component of 
the electric field, ignoring non-linear EM perturbations and invoking the usual gyrokinetic 



2 



ordering for admissible wave forms. This implies that, the Maxwell's equations in his ap- 
proach were actually modified to maintain the ordering properties invoked for the linear EM 
perturbations. These deficiencies were corrected by Beklemishev et at (1999-2005, 0, S I^l) 
who developed a fully covariant theory (derivation of the gyrokinetic transformation) hold- 
ing in principle for arbitrary curved space-times. This theory, in particular, is applicable in 
the presence of non-uniform and non-linear electromagnetic (EM) and gravitational fields, as 
well as arbitrary wave-fields, i.e., non-linear perturbations of the EM and gravitational fields. 
Basic features of their approach, consistent with the allowance of relativistic drifts typical 
of many astrophysical plasmas, were the assumption of non-vanishing parallel electric field 
and the description of particle dynamics in the field-related basis tetrad, allowing a simpler 
formulation of the relativistic equations of motion. In addition, the perturbative theory 
was employed to derive the relativistic Vlasov-Maxwell equations expressed in gyrokinetic 
variables accurate through second order in e (which is the ratio of the gyro-radius/gyro-time 
to the equilibrium gradient length/time), while the wave fields are taken into account in the 
second order in the amplitude. Another remarkable feature of the formulation was to extend 
the customary applicability limits of GKT to include wave fields which can be fast oscil- 
lating in time (oop^fc ~ 1) as well as in space (kp^ ~ 1), without involving the customary 
assumption k\\pi, <C l(feii denoting the parallel wave- vector of the perturbation with respect 
to the local direction of the magnetic field). Nevertheless, the approach of Beklemishev et at 
contained a limitation linked to the case of a small or vanishing background electric field. If 
E and B denote the invariant eigenvalues of the electromagnetic (EM) field tensor F^ V) this 
happens if there results locally E/B <C 1( denoted as weak electric field assumption (WEA); 
see also below) or even it occurs E = 0. For a certain class of problems in plasma-physics 
and astrophysics this limitation may be serious. In particular, this may occur (in relativistic 
plasmas) either if: 

1. The plasma is characterized by a strong turbulence. Turbulence can manifest itself by 
the possible simultaneous presence of perturbations characterizing both the EM field 
and the plasma (to be described either in terms of kinetic or fluid descriptions). In 
such cases, due to strong perturbations, the electric field may locally vanish or become 
much smaller than the magnetic field. 



2. However, difficulties with the perturbative gyrokinetic treatment adopted in 
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0] may occur for slowly varying (or " quasi-equilibrium" ) plasmas and even in the 
absence of wave fields if WEA occurs. This may happen when the standard gyrokinetic 
perturbative theory is carried out at higher orders in the Larmor-radius expansion. 

3. The magnetic field is produced primarily by strong diamagnetic currents. In fact, even 
in the absence of an electric field (i.e., letting E = 0), it is well known that electric 
currents are produced in a magnetized plasma by gradients of the relevant fluid fields 
(in particular, fluid velocity, number density, temperature, magnetic field intensity, 
etc.). 

The physical motivations of the present investigation are related to the occurrence of 
strongly magnetized relativistic plasmas in curved space-time, which may be strongly char- 
acterized by the presence of EM and kinetic turbulence. In the astrophysical context they 
are typically related to the existence of accretion disks, plasma inflows and outflows and 
relativistic jets, all occurring typically close to massive objects, such as neutron stars, black 
holes (BH) and active galactic nuclei (Mohanty [3j). These plasmas are believed to be domi- 
nated by collisionless, fully relativistic dynamics in which the charged particles of the plasma 
interact mutually primarily by means for the very EM mean field which they produce. This 
limitation was overcome by the GKT approach developed by Tessarotto et al. (Ref.l [lo| ) 
extending the earlier formulation 7|,|8|, |9|]. The new approach is based on key features: 1) the 
adoption of a synchronous hybrid Hamilton variational principle to describe single-particle 
covariant dynamics; 2) the introduction of an extended gyrokinetic transformation, which 
includes now also a suitable 4- velocity transformation and is taken of the form 
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)^y i = {r' a ,u' a ), (1) 



where r a and u 13 coincide with the countervariant components of the 4- vectors r (4-position) 
and u (4- velocity), r /a and u' a denote the gyrokinetic variables; 3) the 4- velocity transforma- 
tion is defined in such a way to warrant the fulfillment of the physical realizability condition 
for the 4-velocity, upu 13 = 1, also by the transformed 4-velocity, namely u'^u' 13 = 1; 4) the 
gyrokinetic transformation (JT71) is defined in terms of a power series expansions of the form 



vf 1 = u'>* 



+ J2s s r^r^ + er^(y',e), (2) 

8=1 
OO 

®J2^(y) = u "®^(y',e) (3) 



8=1 



(extended gyrokinetic transformation), where rf (y l ) and vf (y % ) (i = 1,..) denote suitable 
perturbations and the operator © is defined by the 4-velocity addition law [6] which per- 
mits to satisfy the physical realizability conditions for all 4-vectors involved and 
evi(y',e)). This requires that, if and denote two arbitrary 4-velocities, i.e., two unit 
time-like 4-vectors, such that a M a M = 1 and 5 M 5 M = 1, the © operator is defined so that 



defines a 4-velocity too. In addition, here the perturbations rf (y l ) and vf (y l ) are defined 
in such a way that both have - by assumption - a pure oscillatory with respect to the 
gyrophase 0'. In Ref.l this approach was adopted to treat finite Larmor-radius effects on 
particle dynamics in the presence of non-uniform electromagnetic (EM) and gravitational 
fields. Nevertheless, an important issue concerns the extension of the theory (developed in 
Ref.l) to take into account also the possible presence of wave- fields, in principle both EM 
and gravitational. In this paper we intend to extend the applicability range of the GKT of 
Ref.l to such a case. On the other hand since the theory relies on a perturbative expansion 
to determine relevant equations of motion, critical issues are its asymptotic convergence and 
accuracy. In this paper we intend to develop a perturbative solution scheme for the covariant 
gyrokinetic equations of motion for a particle immersed in a strong background magnetic 
field, which exhibits the correct asymptotic convergence properties even in the presence of 
EM perturbations. In detail, the main goals of the investigation are as follows: 

1. generalization of the approach developed in Ref.l to include the treatment of wave- 
fields (EM and gravitational) perturbations based on the introduction of a suitable 
gyrokinetic 4-velocity transformation; 

2. explicit construction of the relevant extended gyrokinetic transformation appropriate 
for the gyrokinetic treatment of linear wave perturbations. In particular, we intend to 
show that, for prescribed wave- fields, the 4-velocity transformation can be reduced to 
the construction of a unique sequence of Lorentz transformations to be constructed in 
such a way to assure a) the fulfillment of the physical realizability condition u^u 11 = 1 
for the 4-velocity; b) the proper convergence of the perturbative expansion also in the 
case of validity of the weak electric field assumption. 



a u © 6, 



a M + b 



(4) 
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These features appear particularly relevant for the treatment of turbulent plasmas in 
the astrophysical context, thus providing an effective tool for use in plasma-physics-related 
problems in astrophysics. 



II. FORMULATION OF CO VARIANT GKT 

Let us briefly summarize the basic steps of the covariant GKT approach based on the 
extended gyrokinetic transformation (G]),©- For this purpose, we adopt here a direct per- 



il!. As 



turbative construction method rather than a non-canonical Lie-transform method 
a starting point, this involves the introduction of a suitable perturbative expansion also for 
the EM 4— potential and the metric tensor g^ v . In particular let us assume that they are 
of the form 

A» = -A"(r a ,£) + Xa^(r a /\), (5) 
g, v = G^(r a 1 e) + \gj lv (r a /\) } (6) 

where e and A are suitable dimensionless real parameters, both to be assumed infinites- 
imal and with £ C A. Here the notation is standard (see Ref.l). Thus, and g^ v 
are respectively the counter-variant and covariant components of the EM 4-potential of 
the metric tensor, while iA M (r a ,e), G M y(r a ,e) and {r u / X) , Xg^ir" / X) denote respec- 
tively the equilibrium terms and the so-called wave-fields, to be identified with suitable, 
rapidly-varying perturbations. All functions here considered are generally assumed to 
be represented by smooth, i.e., in functions. In particular, invoking the gyroki- 

netic transformation (T5]) for the 4- vector r a , it follows that "equilibrium" fields can be 
expanded in Taylor series near the guiding-center 4-position r' a . This implies that, neglect- 
ing corrections of order o(e), ^A^(r a ) and G^r", e) can be approximated respectively as 
^j- = ^j- + Aq [1 + o(e)] and G^ v = G' + sG 1 [1 + o(e)] . Here the primed quantities are 
all evaluated, as usual, at the guiding-center position = A M [r' a ) , G' = g^ {r' a ) , 
Aq = r^d' a A'^ and G l v = r^d^G' ), while the perturbations (of relevant fields) depend 
generally also on the gyrophase angle 0' through the perturbation erf (y l ) . The perturba- 
tive scheme converges, at least in an asymptotic sense, only if one assumes the validity of 
a suitable ordering conditions (Larmor radius ordering), implying that EM contributions 
must dominate, in a suitable sense, over inertial ones in the variational functional. For- 
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mally this corresponds to replace the electron electric charge e by e/e, leaving the rest 
mass mo unchanged and considering the 4- velocity of order o(e°) . To obtain the GKT for 
particle dynamics the relevant variational functional S(y) must be expressed in terms of 
suitable gyrokinetic variables, y l , which include in particular the gyrophase <fr describing 
the fast gyration motion of particles around magnetic flux lines. As shown in Ref.l a 
convenient formulation can be achieved by identifying the variational functional with the 
action S(y) = J 2 7(y) + dF, where y denotes the set of variables y = (r M , u^, A) and 7 is the 
differential 1-form 7(y) = g^ v [qA u + u v \ dr^ + £ [g f _ lu u tJl u u — 1] , with £ a Lagrange multiplier 
and q = Ze/mo the normalized electric charge. It is immediate to show that phase-space 
trajectory of the point particle with 4-position r u and 4-velocity u v results, by construc- 
tion, an extremal curve determined by the synchronous variational principle SS(y) = 0, 



where the synchronous variations 5r^,8u^,8^ are all considered independent [lOj. Indeed, 
the physical realizability condition u u u v = 1 is satisfied only by the extremal curve. To 
construct the transformations ((2j)-^3^ corresponding to (JHJ)- flHj) the perturbations rf (y l ) and 
vf (y l ) must be suitably determined. As usual, u' a can be expressed in terms of suitable 
independent (generally hybrid) gyrokinetic variables, one of which includes by definition the 
gyrophase 0'. This is obtained by projecting u' a along the four independent directions defin- 
ing the so-called fundamental tetrad unit 4- vectors (with a = 0, 3), hereon also denoted 
(Z /At , Z //At , Z M , r^) [or in symbolic notation (1, 1', 1", r)], where Z M and are respectively 

space-like and time-like orthogonal unit vectors. In particular, the latter can always be 
identified with the EM fundamental tetrad, i.e., the basis vectors of the EM field tensor, 
in this case to be evaluated at the gyrocenter position, F' = F^ u (x' a ), with eigenvalues B 
and E. In particular, the space-like vectors e£ , e^ 3 ) are assumed to satisfy the eigenvalue 
equations F'e% = Be^, F' e% = —Be2 il . Here e aAt = g'„ v e v a are the covariant components of 
and the labels a assume a tensor meaning in the tangent space [so that in terms of the 
Minkowskian tensor 77 = diag (1,-1,— 1,-1) there results e a = rj ab eb^\- 



III. CONSTRUCTION OF THE EXTENDED GYROKINETIC TRANSFORMA- 
TION 

To illustrate the key new features of the present approach let us point out the role of the 
extended gyrokinetic transformations (T5]) and ([3]). GKT is a perturbative theory which can 
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in principle be carried out at any order - both in e, and A - provided the perturbations (i.e., 
eff and evf ) are suitably small, i.e., in particular there results uniformly, in an appropriate 
subset of phase-space, 

eff ~ o(e). (7) 

This requirement is clearly necessary for the convergence of the gyrokinetic perturbative 
expansion. Hence its must be viewed as a consistency condition for the validity of GKT. To 
elucidate the point, it is sufficient to consider the perturbative formulation of GKT based 
on an expansion in A which is accurate only to leading order in the dimensionless parameter 
A. As indicated above, A represents the magnitude of the wave fields defined by Eqs.flS]) 
and (jSJ). Hence in this case a linear approximation in the wave-fields is actually considered 
for the perturbative expansion. If the velocity perturbations in Eq.Q are neglected, one 
can prove that to assure that the differential 1-form 7 results gyro-phase independent, the 
covariant formulation of GKT requires that to leading order in A the following two constraints 

Q 

equations must be satisfied [7J]: 

«;+ W„ - qr\F^ = (8) 

>4 + ^-^)§£}~ = fP. <») 

where R represents in principle an arbitrary gauge function (which can always be set R = 0). 
Assuming that EB 7^ and ignoring possible 4- velocity corrections [see Eq.flSJ)], Eq.(jHJ) 
delivers a formal solution for the perturbation r" (Larmor-radius 4- vector) which to leading- 
order in A yields r\ = — cos <p' + Z£ sin <p'] +XD h, i*a li . Here D' Ufl is the inverse Faraday 
tensor which reads D'^ = F'~ l = \(l v r» - 1^t u ) - ±{l' v l"v - l'H" v ). This expression 
manifestly diverges in the limit E — > 0. In particular, for \E\ sufficiently small it is obvious 
that the second term on the r.h.s. can become arbitrarily large thus violating the ordering 
condition (J7j) (weak electric field assumption) . To assure the proper convergence of GKT also 
in this case it is manifest that an extended formulation of GKT is required. In particular, 
as indicated in Ref.l, it is actually necessary to introduce a perturbative expansion also for 
the 4-velocity. This can be conveniently represented in terms of the 4-velocity perturbative 
expansion (j3J) and by imposing the 4-velocity addition law (j3J). In such a case the constraints 
flU) and Qj are replaced by the three equations: 

K © KJ~ + X K - <i F ^ = °> ( 10 ) 



{ ( W © + A < - = 0> (n) 

{[u' © Am' 01 ] 2 - [wj © Au[|J 2 - [«;' © Xw[] 2 - l}~ = 0. (12) 

Its is immediate to prove that - to leading-order in A - these equations can be satisfied 
identically for arbitrary values of E, by suitable definition of the 4- velocity perturbation 
Xvip, defined in such away to realize appropriate velocity cancellations for the wave-field 
perturbations Xqa^ occurring in the previous equations. In particular it can be proven that 
this result can be achieved by means of the following two Lorentz transformations: 

• First transformation (Lorentz boost, space-time rotation): Let us first introduce a 
space-time 4-rotation acting along the directions (r, 1) and let us impose that the 
spatial component of the transformed velocity cancels identically the corresponding 
component of the wave perturbation a. The 4-rotation takes the form: 

A 

ujl = («j| © Ai>i||) cosh a + (u' © Xv' 10 ) sinh a 
a (13) 
u'q = © Ai^||) sinh a + (u' © Xv' 10 ) cosh a 

A A 

where uj| e u' represent the two surviving components of the 4-velocity after the 
transformation. 

• Second transformation (Lorentz- boost, space-time rotation): Let us introduce now a 
second space-time rotation acting along the two remaining directions (l',l",r), and 

A 

which leave unchanged the parallel component u'» (along the unit 4-vector r) already 



determined by the first Lorentz boost. We therefore define the second 4-rotation: 

(14) 



_J A 

w — (w' © Xw[(y)) cosh {3 + u' sinh (3 

— A 

u' Q = u' cosh (5 + (w' © Xw[(y)) sinh (3 



After the second transformation the 4-velocity takes the (desired) final form 

m m = w' % cos(0') + l'l sin(0')] + u\l^ + u'^. (15) 
As a result one finds that the correct definition of the Larmor 4-vector is 

r\ = -D'^ + XD'^a'^ = (16) 



-2^" [ft cos 0' + Zjffln0O] + A^'^ (a^ + a 3 Zj) , 



while the gyrokinetic differential 1-form (which is by construction gyro-phase independent) 
reads 

7(y') = (f 4 + < © X( i%) dr " + \ (K © A A t|) + 

(18) 

This expression is formally similar to that determined in Ref.l, except for the explicit ap- 
pearance of terms containing the 4- velocity addition law (j4]) on the r.h.s. of Eq. (fl7l) . which 
have been introduced to preserve their proper physical interpretation (as 4-velocities) in 
the corresponding Euler-Lagrange equations. This implies that non-linear corrections in the 
wave-fields are, actually, taken into account in the gyrokinetic differential 1-form. This result 
shows the (formal) intrinsic simplicity of the perturbative approach to GKT here developed. 



IV. CONCLUSIONS 



In this paper we have formulated a covariant gyrokinetic approach for single-particle 
dynamics utilizing the variational approach developed in Ref.l. We have shown that: 

• the theory goes beyond the domain of validity of the treatment previously considered 
3, [§], E3] and includes the case in which the electric field can locally vanish or result 
much smaller than the magnetic field; 

• the treatment includes the effect produced by the presence of wave-fields (produced 
by EM and gravitational perturbations). In particular, we have shown that this gen- 
erally requires the introduction of an extended gyrokinetic transformation based on 
the introduction of an appropriate 4- velocity perturbative expansion; 

• a relativistic addition law has been invoked for the 4-velocity perturbative expansion; 

• non-linear corrections in the wave fields have been take into account in the pertur- 
bative theory to assure that the condition of physical realizability is satisfied by the 
transformed (gyrokinetic) 4-velocity. 

The simplicity of the approach and its applicability to a wide range of possible physical 
situations make the present theory particularly suitable for applications to the investigation 
of relativistic plasmas, potentially both in astrophysics and laboratory plasmas. 
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